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INTRODUCTION 


The effect of temperature in acoustics is usually analyzed by allowing the speed of 
sound to be a function of spatial position 1 ’ 2 . This type of analysis ignores the explicit 
effect of the temperature gradient. It is known that the acoustic perturbation equations 
contain explicit derivatives of the thermodynamic state variables, such as pressure and 
density, but these derivatives may often be ignored because they are small on the scale of 
an acoustic wavelength. Explicit gradient effects have been included in studies of sound 
propagation in the atmosphere 3 ’ 4 where long wavelenghs are of interest. These wavelengths 
are called infrasonic, where the temporal frequencies are on the order of 1 Hz and less. 
The temperature gradient in the atmosphere is typically on the order of several kelvins 
per kilometer so that the wavelength must be long in order for the temperature to vary 
significantly over a wavelength. Gradient effects may be neglected in studies of audio- 
frequency (20Hz-20kHz) sound in the atmosphere. 

The purpose of this paper to is analyze problems where the temperature gradient is 
strong and cannot be neglected. One problem considered is a simple rectangular cavity 
with the temperature being a function of only one spatial coordinate, the one representing 
the height of the cavity. This situation may occur, for example, within the structures of 
high-speed aircraft where the skin is heated by a supersonic or hypersonic boundary layer. 
The cavity could be the interior of a Thermal Protection System (TPS) 5 structural panel 
such as shown in Fig. 1. The outer skin can be up to 1350 kelvins while the inner skin 
should be relatively cool, say 300 to 400 kelvins. The temperature gradient within such 
a panel would then be 10 to 20 kelvins/cm, or six orders of magnitude larger than the 
typical atmospheric gradient. More importantly, it is large for wavelengths of interest in 
the analysis of the structural panels. The resonant frequencies of the cavity are important 
parameters in the dynamic response of a panel over the cavity. If the temperature gradient 
significantly affects these resonant frequencies, then the gradient is important in predicting 
the panel response. This paper will define the effect of a temperature gradient on the 
resonant frequencies of the cavity. 

The panel shown in Fig. 1 is filled with insulating blankets. The acoustic properties 
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of such blankets may be described by the Attenborough 6 model in the case of moderate 
temperatures, but the inclusion of these properties would complicate the present analysis 
by introducing the fibrous materials’ parameters. Also, the materials would reduce the 
resonant effects studied here by way of increasing the damping of the acoustic waves. 
A simpler model, shown in Fig. 2, is used here in order to focus on a single effect — 
the temperature gradient. The cavity is assumed to have perfectly rigid walls with the 
temperature varying linearly between the bottom and the top. A perfect gas at constant 
pressure fills the cavity so that the density varies inversely with temperature. 

The essential features of the resonant cavity response may be verified by an experiment 
on a waveguide. The waveguide, shown in Fig. 3, is like the cavity except that the width 
is equal to or less than the height. The restricted width eliminates cross modes which 
depend on the spatial dimension y and thus simplifies the analysis by restricting it to the 
x-z plane. 

A second problem where the strong gradient analysis has application is the response 
of a recessed pressure sensor. Sensors are installed in recesses when the wall temperature 
in high-speed boundary layer flows are above the operational limits of the pressure sensors. 
When this technique is used, a transfer function must be developed to relate the sensed 
pressure fluctuations to the wall pressure fluctuations. Tijdeman 7 ’ 8 , has given the transfer 
functions in recesses where the temperature is constant. These transfer functions show 
resonant peaks which depend on the cavity depth. The analysis here may be used to show 
the effect of a strong temperature gradient on the transfer function for a recessed sensor. 

The first section of the paper is devoted to developing the acoustic wave equation 
with a temperature gradient effect. More general equations have been developed than the 
one used here. Pridmore-Brown 3 has given the perturbation equations including pressure, 
density, and velocity gradients. Grossard and Hooke 4 give completely general equations in 
a rotating coordinate system to represent the earth's atmosphere. The specialization here is 
intended to emphasize a single effect. Consequently, the pressure is assumed to be constant 
and there is no flow. These simplifications produce a wave equation which is similar to the 
Helmholtz equation used in elementary acoustics. The equation is generalized, however, 
by the inclusion of the temperature gradient in a coefficient. This generalized equation, 
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called the strong gradient equation, shows why the resonant frequencies of a cavity will be 
altered by the gradient. 

The second section of the paper is devoted to the response of a simple rectangular 
cavity to a monopole source. A set of orthonormal basis functions is constructed such 
that the gradient-dependent boundary conditions are satisfied. Asymptotic formulas are 
derived for the effect of the temperature gradients on the modal wave numbers and phase. 
The Galerkin 10 method is used to express the cavity response in terms of these modes. 
All integrals involved in the Galerkin procedure are evaluated by exact formulas 1 . Second- 
order accurate asymptotic formulas are developed for the natural frequencies, modes, and 
forced response of the cavity. 

An experiment is suggested in the third section to verify the effects predicted here. 
The experiment would use a rectangular waveguide which is simpler than the rectangular 
cavity. The width of the waveguide is selected to be small in comparison to the height so 
that fewer modes will propagate. The length of the waveguide is large in comparison to the 
height. Standing waves are excited by sources at one end of the waveguide and monitored 
by microphones at the other end. Temperature and temperature gradient are controlled 
by heated and cooled conducting plates forming the top and bottom walls, respectively, of 
the waveguide. The other sides are near-perfect insulators so that temperature variation 
is one-dimensional in accordance with the theory. Excitation frequency is varied with 
excitation voltage held constant so that the microphone output represents the frequency 
response function of the waveguide. Comparison of the response functions for different 
combinations of plate temperature shows the temperature gradient effects. Predictions of 
the resonant peaks are made to simulate the proposed experiment. 
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SYMBOLS 


ai 

A n > 

B u‘ 

c« 

c 

Crp 

C P 

f 

9 

H 

k 

k 

L 

P 

Q 

T 

T 

t 

u 

W 

x,y,z 

Greek 

a 

0 

1 

6 

€ 

e 


mode normalization constant 
array of weak gradient integrals 
array of strong gradient integrals 
Array of constants 
sound speed, m/s 

sound speed at average temperature, m/s 
specific heat at constant pressure, J / Kg ■ K 
frequency, Hz 
mode function 

height of rectangular cavity, m 

unit vector in z-direction 

wave number, m~ l 

length of rectangular cavity, m 

pressure, Pa 

source strength, s -1 

temperature, K 

average temperature, K 

time, s 

velocity, m/s 

width of rectangular cavity, m 
cartesian coordinates 


frequency-shift constant 

acoustic admittance 

ratio of specific heats 

Kronecker delta fun< (ion 

dimensionless gradient parameter, AT/T 

phase angle, rad 
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K 

£ 

p 

T 

$ 

U) 

Subscripts: 

l 

m 

n 

0 

P 

T 

x 

y 

z 

T 

Special Symbols: 

/ 

A 

/ 

/ 

/' 


thermal conductivity, W/m ■ K 
dimensionless coordinate, z/H 
density, kg /in 3 
thermal gradient 

modified acoustic potential f 
acoustic potential 
circular frequency, rad/s 


mode shape in z direction 
mode shape in y direction 
mode shape in x direction 
origin 

constant pressure 
constant temperature 
x component 
y component 
z component 
thermal gradient 

acoustic quantity / in time domain 
acoustic quantity / in frequency domain 
average /, or based on the average / 
derivative of / with respect to it’s argument 
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STRONG TEMPERATURE GRADIENT EFFECTS 
Acoustic Equations 

The following analysis is based on the essentially-standard equations of fluid mechan- 
ics, namely the continuity, momentum and energy equations. The continuity equation is 
taken to be 

( 1 ) 


§7 + » v ' u = pQ 


The source term on the right of the continuity equation represents a fluctuating volumetric 
source distribution of strength Q, which is small in some sense, representing the magnitude 
of all acoustic terms. The momentum equation is written without viscous effects included. 

DU Vp 


Dt + p 


0 


(2) 


The energy equation is written with heat conduction effects, but without viscous effects. 
The heat transfer term will be neglected in the acoustic equations on the grounds that the 
thermal diffusivity is small relative to a typical wqavelength and period of the phenomenon. 
The energy equation is 

D n n D /i 

( 3 ) 


^ =1 ?^ + h _l )KV »T 
Dt p Dt v ' ' 


In a system without flow, the velocity is decomposed into a null vector and an acoustic 
vector. 

U <- 0 + u (4) 


The steady momentum equation then requires the pressure gradient to be zero. 

Vp 

— = 0 


( 5 ) 


The steady energy equation requires that the temperature satisfy the Laplacian equation. 


V 2 T - 0 


( 6 ) 


The density is decomposed into steady and fluctunl jne parts to form the acoustic continuity 
equation. 

p - p + p (7) 
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The acoustic part of the continuity equation becomes 


dp _ _ ~ 

-- + u-Vp + pV-u = pQ 
dt 


Similarly, the acoustic part of the energy equation is 


( 8 ) 


dp (dp \ 

p Tt =lp U + '’J 


(») 


where the heat transfer term kVT has been neglected. The energy and continuity equations 
combined relate the fluctuating pressure to the fluctuating velocity and the acoustic source. 


~ 1 dp ~ 

V-ii H — ^ = Q 

7 p dt 


( 10 ) 


If an acoustic potential is defined by 


P = P 


d* 

dt 


( 11 ) 


then the continuity-energy (10) equation becomes 

„ _ Id 2 * ~ 

The acoustic part of momentum equation is 

<?u Vp 

+ = 0 

dt p 

and the acoustic velocity in terms of the potential as 


d 2 u d 2 * VTd 2 * 

= -V^r + 


( 12 ) 


( 13 ) 


dt 2 


dt 2 T dt 2 


( 14 ) 


where the gas law, p ~ pRT , has been used to eliminate density in favor of tempera- 
ture. The acoustic wave equation is formed by solving the momentum equation for u and 
substituting this solution into the continuity-energy equation (12). 

Harmonic solutions of the form 


u - ue~ lwt 


( 15 ) 


7 




(16) 


convert the continuity and momentum equations into 

V-u — k 2 4> = Q 


(17) 


and 


„ VT » 
u = -V* + — $ 


(18) 


Wave Equation 

When the vector VT is null, equations (17) and (18) combine to form the elementary 
acoustic wave equation 

V 2 tf + A; 2 * = -Q (19) 


When the temperature is variable, the momentum equation (18) gives an additional 
acoustic velocity component in the direction of the temperature gradient. The continuity 
equation then gives a wave equation with the temperature gradient as a variable coefficient. 
The heat conduction condition, equation (6), allows the wave equation to be expressed 
without V 2 T. 


0 A VI A 

V 2 * - — Vtf + 




( 20 ) 


The gradient of the acoustic potential is eliminated by introducing a modified potential $ 
with an integrating factor. 

* - V'7’/T 0 $ (21) 


The resulting wave equation then depends on the magnitude of the logarithmic gradient 


of the temperature. 


V 2 $ + 




( 22 ) 


The gradient term in the wave equation can be expressed as a frequency parameter. The 
product of the speed of sound and the logarithmic derivative of any variable has the 
dimensions of frequency so that the frequency pa m motor is defined as 


c 

L ° r = 2 


VT 

T 


= |Vc| 


(23) 
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and will be called the temperature gradient frequency. The effective frequency in the 
wave equation is the hypotenuse of a frequency triangle formed with the actual excitation 
frequency u and the temperature gradient frequency ui T . The gradient frequency makes 
the strong gradient equation look very much like the weak gradient equation 


V 2 $ + 



(24) 


Boundary Conditions 

The boundary condition for a locally-reacting surface is expressed as a proportion 
between the acoustic pressure and the acoustic velocity normal to the surface. 


P = Zu n 


(25) 


The acoustic pressure is proportional to the potential function. 


p = — itjjpAQ 


(26) 


The normal velocity depends on the normal gradients of both the temperature and the 
potential. From equations (14) and (20) 


u r 


To i dn 2 T dn 


(27) 


The boundary condition for the potential contains the effect of the temperature gradient 
normal to the wall. 


_ /J_ 
~ \2T 


dT 


4- ik/3 1 $ 


(28) 


2 T dn 

The effect of the temperature gradient normal to the wall is analogous to an acoustic 
admittance (3 T . 


0T — 


— i dT -i. dc 


2 k.T dn te dn 

which is purely imaginary and thus introduces no dissipation. 


(29) 
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ACOUSTIC RESPONSE OF A RECTANGULAR CAVITY 


Governing Equations 

The strong-gradient wave equation will be used to predict the response of a rectangular 
cavity to a monopole source of strength Q s . Place the origin of coordinates in the center 
of the cavity as shown in Fig. 2. With the temperature gradient in the positive z-direction, 
the temperature is given by 

T = T (1 + tz) (30) 


The thermal gradient frequency is a function of the gradient r. 

M 


U r = 


C 


2 y/l + TZ 

A 

The acoustic potential $ is the solution to the following wave equation 

A 

-Qs 


(31) 


+ 


^ 2 /4 + 

1 + TZ (1 + TZ)‘ 


$ = 


6(x - x 3 )6[y - y 3 )S(z - z 3 ) (32) 


y/T+Tza 

The solution will be developed for the case of rigid walls such that the acoustic velocity 
normal to the wall is zero. The velocity is a linear combination of the potential gradient 
and the temperature gradient vectors. 


u = —y/l + rzV$ + k ( - = ) $ 

\2 vT + tz ) 


(33) 


The boundary conditions for the acoustic potential in the rectangular cavity are homoge- 
neous and depend on the gradient in the corresponding direction. 


dx 


= 0, x 


' 2 

W 

- 0 , y = ± — 

dy 2 

, . t - H 

( l + rz hi = i*' z = ± T 


(34a) 

(346) 

(34c) 


The wave equation (32) and boundary conditions (34) define 4> within the rectangular 
cavity. The x-boundary conditions (34a) ar** satisfied by a s*4 of modes <i„(k rr ,x) defined 

below. 

~ \/2,cos(k xn x 0 n ) , ti — 0, 1 , 2, . . . (35a) 
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. U'K 

kxn ~T 

(356) 

II 

O 

(35c) 

. JO, if n is even; 

n 1 7 t/ 2, if n is odd. 

(35d) 


The integral of g 2 (k xn x) over the interval (— L/2, L/2) is equal to the length of the interval. 



x n'%)dx — <5nn' L 


(36) 


The potential function is expanded in a double series of modes whose amplitudes depend 
on the vertical position z. 


OO OO 


$(x,y,z;w) = EE ( 2 i k-Offn (&xn ;c )<7m(^j/ml/) 


n=0 m= 0 

Each function $ nm (z;u>) is defined by an ordinary differential equation. 


d 2 $ r 


dz 2 


+ 


«V4 , T U± 


1 + TZ (1 + Tz) 


2 "nm 




Ze) 

y/l + TZ 


The horizontal mode numbers are 


(37) 


(38) 


JL/ _ L* I L>* 

A/ nm ~ a xn ' a ym 


(39) 


The generalized source depends on the monopole source strength, and the values of the 
mode functions at the source position. 

Qnm ” jjyy 9n{kxn% a)ffrn(fcym!/ 3 ) (40) 

Finally, the boundary condition (34c) carries over to each mode in the same form. 

/-i \ d^ n77 i t - H 

{1 + rz) -dT^ (4)| 

A rough estimate of the natural frequency of a mode in a thin cavity is given by 
assuming 

£ + 7 = k nm (42) 

T ^ 
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This simple estimate indicates that the natural frequency will be less than the natural 
frequency w nm = Cfk nm of a cavity without a temperature gradient. 


Wnm — 




(43) 


The effect of the gradient-dependent boundary conditions on the acoustic modes is devel- 
oped in the following section. 

Basis Functions 

Basis functions in the vertical dimension are generalizations of the cosine functions in 
the horizontal dimensions. 

Ql(kziz) = aicos(k zl z - 0 zl ) (44) 

These functions must satisfy boundary conditions which depend on the temperature gra- 
dient. 

r H 

(1 + Tz)g' l (k z iz - 9 Z i) = -gi(k z iz - 0 Z j), a = ±— ( 45 ) 

The boundary conditions yield the following characteristic equation for the wave number 
k z i . 

tan k zl H 2e 2 ^ 


k zl H (4 - c 2 )* 2 ,# 2 + e 2 

e = tH — AT /T (47) 

The series for tan x/x contains only even powers so that ( k z iH ) 2 is a function of (AT /T) 2 . 
The branch is selected such that the two variables have the same sign when they are real- 
valued. The series also provide asymptotic formulas for the wave numbers. 


k z0 H 


k zi H 


=in { i+ d 


e 2 + 


e + 


) 


(48a) 


(486) 


The lowest order wave number is proportional to the temperature gradient, but the higher- 
order wave numbers are functions of the square of the temperature gradient. 

The phase of the modes is found from eider boundary ">nditi»n. 

(2 ± t)k z iHsm k z [H/2 :t ecos k z iH /2 


tan 9 z i — ± 


(2 ± e)k z \Hcos kJH/2 =F esin k z iH / 2 


(49) 
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The dominant term in the asymptotic formula for the phase of the modes shows that the 
modes approach the simpler sine and cosine functions of the zero-gradient case. 


„ 7 T 1 , 

0 ? o - — T — c T 

4 24 


= f 0, l even \ e 
zl \ 7r/2, l odd j 2ln 


(50a) 


(506) 


Thus the effect of the temperature gradient is to give the fundamental mode a small wave 
number which is proportional to the gradient. The fundamental mode becomes slightly 
nonuniform, varying slowly with distance, because k z o is proportional to c. The higher 
modes have wave numbers shifted by the square of the gradient, but the phase of the mode 
is shifted in proportion to the gradient. The principal effect is to translate the mode in the 
positive or negative direction, depending on the sign of the gradient, as shown by equation 
(50b). The modes are neither even nor odd functions, but have a mixed character because 
of the gradient. 

The normalization constant a/ depends on the wave number and phase of the mode. 

- 1/2 


ai 


= v^( 


sin k z iH 

1 + ~ ""£ — r T~ co * 20 i 

k z iH 


) 


(51) 


The expansion of equation (51) shows that the effect of the gradient on the normalization 
constant is of order e 2 . 


a/ 


o = \/2 ( 

= v/2(l 


1 + s ( 


•) 


4/ 2 7T 2 


(52a) 

(526) 


Galerkin Solution 

An approximate solution to the modal wave equations (38) and (41) can be made by 
extending the mode expansion to the vertical direction and selecting the mode amplitudes 
by the Galerkin criterion 9 . The result is a set of algebraic equations which are coupled in 
the vertical index but uncoupled in the horizontal indices. 


E 


w 2 , r 2 „ 

-~ 2 -All' + -7B1V 
4 

i '= 0 L C T 4 


^nml' ( ^nm ^ zl ) ^ n ml — Qnml 


(53) 
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(54) 


a - JL [ H ^ 2 srtkf&hw {kzt ,z ) 

W ~ H]_ h/2 (1 +TZ) 

R 1 f H/ 2 gi{kziz)gv{k z i>z) 

B "' = hL 


dz 


dz 


(55) 


-H/ 2 l 1 + rz ) 2 

The integrals in equations (54) and (55) are expressible as sine and cosine integrals. 
Alternatively, power series may be used to develop approximations for the effect of the 
temperature gradient. Second-order- accurate approximations will be used here to clarify 
the effect of the temperature gradient on the modal functions and cavity response. Since 
Bu> is multiplied by r 2 /4 in equation (53), only the zero-order approximation, Sw is needed 
for these integrals. The second-order response equation is thus 


(^Iml + e2 ^ 2 U/U ^) $/ - W 2 ^ Au'$l' = c 2 f Qi 


where 


ui 


-{t 


/ = o 

l > 0 


(56) 


(57) 


and Aui is approximated by the first three terms of a power series in e 

1 £ 2 II 

Au> = 6ir + (.Ant + — A^/ 


(58) 


The frequencies u> nm i are the resonant freqencies of a cavity with constant temperature 


- 2 

^nml 



tin 

W* n = 

Cf - 


mu 

tiym — 



In 

Qlz = 

Cf H 

- Q ln 

+ Q l>n 


-2 


(59a) 

(596) 

(59c) 

(59d) 


It is understood that oo z = u> z i. In the following *'qnal ions, the subscript (n.m) will be 
suppressed for brevity. When Qi appears, it. means u’ nm |. Equivalent interpretations are 

A A 

to be used for the symbols and Qi. 
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The modes and their integrals will be given terms of a dimensionless coordinate £ = 
2 / H which has the range (-1/2, 1/2). The second-order approximations for the modes are 

go(kzoH0*l + ^-£j (60a) 

sin/jr^ ,1 = 1,3,5,... 

( 606 ) 

gi ( k z tH £ ) « \/2 cos In £ + cosing - e 2 ^^£sin/7r£ + cos In ,1 = 2,4,6... 

(60c) 

Equations (60) shows the effect of the temperature gradient on the acoustic modes. The 
modes are no longer even and odd functions in correspondence with their index. Each 
mode has an odd (in £) part, which is proportional to the temperature gradient, and each 
even-numbered mode has an odd part. The appearance of the functions £sin(ln£) and 
£cos (/tt£) in the second-order gradient effects indicates that a Fourier series approximation 
to the strong-gradient solution would be accurate to only first order in gradient. 

The first derivative A' u , are found to be zero unless l + 1' is an odd number. 


/ 2 

9l {k z iH£) « \/2sin/7r£ — e— — cos/7r£ + e 2 

2t 7T 


V2 

2ln 


£ cos In ^ 


3y/2 

8n 2 l 2 



ill,, =0, l + l' = 0,2,4... 

(61a) 


A',o=(-l) <,+,)/2 ^, ( = 1,3.5... 

(616) 

A IV — (' 

+ "* 0 ^ ) il + l =1,3,5... 

’ n 2 (/2 _ j'2) 2 ’ \ / > 0, l' > 0 

(61c) 

The second derivatives are zero unless l + /' is an even number 



1 » 

2^oo — 0 

(61d) 


2^ = <- 1),/2 2(V- ' = 2 ' 4 ' 6 "' 

(61e) 


1 „ ( l 1 \ 

-An, = ( — + 1 = 1.2.3... 

2 11 \ 12 l 2 n 2 I 

(61/) 

1 A " 
2 Aw ~ 

( Dtl+i ')/ 2 4 V +/2 ) f/f/' =2,4,6... 

1 ' n 2 (/ 2 _ /' 2 ) 2 ’ \l > 0 , /' > 0 

(61?) 
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Natural Frequencies and Modes 


Homogeneous solutions of equations (56) are possible at resonant frequencies uji which 
depend on the temperature gradient parameter e. These resonant frequencies are solutions 
of the infinite determinant equation 

|u; 2 j4u' — f>w (u» 2 + e 2 s/tt> 2 )| = 0 (62) 

When this determinant is expanded, it is found that the cofactor of any element of order 
e is also of order e. Consequently, the resonant frequencies are functions of e 2 and can be 
approximated as 

<Ajf « uf (l + e 2 a;) (63) 

The corresponding natural modes, designated by xu are given as 

xu = x{ t 0) + ex\p + e 2 xj 2) (64) 

Since the homogeneous solution (64) may be multiplied by an arbitrary constant, there 
is no loss of generality in taking x^ = 1 , = 0 and xj^ = 0. The remaining elements 

of the vectors are found by substituting the assumed solution (64) with the series (58) for 
Aw into the homogeneous equation from (56) and collecting coefficients of each power of 
e. The zero order vectors are the Kronecker delta function 

x\ft = Su> (65) 


The first order vectors are proportional to the first derivatives Au> 


x (1) - D 2 
X ll' ~ U l 


°, , 
A u> 
l ' 3 -/ 3 


l =/ 

l^l' 


where P{ is the reduced frequency for the mode 


P 2 = 


oh 


nH 

17 


( mH \ 2 fJ 

+ IF +, ‘ 


( 66 ) 


(67) 


The second-order vectors are found to depend on the second derivative A ' a , and on products 
of the first derivatives. 


x 


( 2 ) 

IV 



1 -4 . p 

2 P-J' 2 V* 


0, 2 C„. 


i?i 


( 68 ) 
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where 


'll’ 


A 00 

-(§) E 




4 * (/' 5 

3*1* 


P) 


( 69 ) 


The frequency-shift constant is found to be a function of pf from the second order equation. 

a, = (£j CuPf - ±A U , + Pi > 0 (70) 

Returning the subscripts n and m, the final equation is given for the natural frequencies of 
a rectangular cavity with a strong temperature gradient. These equations are the principle 
result of this paper 



w nml & nml 

(AT 

1 + <*nmi I ^ 

)’] 

(71«) 


&nmO — | 

^2\ 4 r *l m0 
W 00 ^001 


(716) 

&nml — ^ 

2 )Vr-( 

* J “OQI V 

' 1 1 \ 
,12 + l 2 7T 2 ) + 

3 / > 0 
4 **0* ml ’ 

( 71c) 


The constant Coo is negative, roughly -1/2. Consequently, the effect of the temperature 
gradient is to shift all resonant frequencies u> mn o to lower values. The change in the square 
of the frequency is proportional to the square of the temperature gradient. For all higher 
modes, l > 0, the resonant frequencies may shift in either direction. The shift of frequency- 
squared is again proportional to gradient-squared, but the constant of proportionality is 
a quadratic function of the square of the reduced frequency for the mode. This reduced 
frequency is equal to or greater than the mode index / and, in the case of equality, it 
can be shown that the quadratic function has a value less than zero so that resonant 
frequencies shift to lower values. The quadratic function is greater than zero for large 
reduced frequencies because the constant Cu is positive. There is a value of the reduced 
frequency Pf , the greater root of the quadratic, which is greater than the index l. Each 
mode whose reduced frequency is in the range [/, iq ) will shift to lower resonant frequencies 
with increasing temperature gradient, and each mod* 3 with reduced frequency in the range 
(tq*, 00) will shift to higher frequency. 

ACOUSTIC WAVEGUIDE EXPERIMENT 
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A rectangular channel with a one-dimensional temperature distribution can be con- 
structed as shown in Fig. 3 to investigate the effect of strong temperature gradient on 
acoustic modes. The width is made small compared to the height so that cross modes, 
m > 0, respond only at high frequencies. It is necessary to cool the lower surface while 
heating the upper in order to keep the temperature extremes within the range of available 
acoustic pressure sensors. A set of five temperature combinations is given in Table I which 
gives evenly-spaced values of (AT/f) 2 from zero to 1. The average temperature is held 
fixed at the boiling point of water, 393.15 K. With this choice, the lowest temperature is 
196.6 K, slightly higher than the sublimation point of dry ice (CO 2 ) which is 194.4 K. The 
highest temperature is 589.7 K. The height of the channel is chosen to give a convenient 
reference frequency of 1 kHz. Based on the speed of sound at 393.15 K, which is 397.49 
m/sec, the height is selected to be 19.87 cm. The length is selected to be a nonintegral 
multiple of height, 5/2, to avoid low-frequency duplication of longitudinal and vertical 
mode resonant frequencies. The length is set at 49.69 cm. These dimensions give temper- 
ature gradients up to 20 K/cm, roughly the same size as those estimated for the thermal 
protection system panel in Fig 1. 

Table II shows resonant frequencies up to 4 kHz for this cavity. The corresponding 
frequency-shift constants are shown in Table III. It is clear from Table III that the largest 
effect of the temperature gradient is on the plane wave (n, 0,0) modes and the first antisym- 
metric (n,0,l) modes. Fig. 4 shows two examples of these modes, the (5,0,0) and (5,0,1) 
modes. 

If the acoustic drivers in Fig. 3 are operated in phase, only the plane wave modes will 
be excited directly. The antisymmetric modes will not be excited because the excitation 
is symmetric. The drivers are located at z — ±H/4 which are the node lines of the (n,0,2) 
and higher order symmetric modes. These positions eliminate the higher order symmetric 
modes. Of course, all modes will receive some excitation as the temperature gradient is 
increased because the modes are no longer perfectly symmetric or antisymmetric and the 
node lines will be shifted. The drivers may be operated 180 degrees out of phase to excite 
the antisymmetric modes with only small excitation of the symmetric modes. 

Resonant frequencies are given in Tables IV and V for the plane-wave and first an- 
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tisymmetric modes, respectively. The (10,0,0) mode has imaginary frequencies in the 
extreme case where AT /T = 1, but this is outside the limit where the asymptotic theory 
is valid. In any event, the tendency of increasing gradient should be to decrease the reso- 
nant frequencies of the plane-wave modes so the waveguide response should be as shown 
in Fig. 5 with symmetric excitation. With no temperature gradient, there should be 6 
resonant peaks between 0 and 2.5 kHz. When the temperature gradient is increased to 
AT /f — 1/2, the lower resonant peaks should be unaffected, but there should be a percep- 
tible downshift of the 4th, 5th and 6th peaks. When the temperature gradient is increased 
to \/2/2 there should be a more pronounced shift and a 7th peak should appear. 

Fig. 6 shows the anticipated response of a waveguide to antisymmetric excitation. The 
response is shown from 2 to 4.4 kHz because the resonant frequencies should shift upward 
in this range. Without a temperature gradient, ther are five resonant peaks, the lowest at 
2236 Hz and the highest at 3736 Hz. As the temperature gradient is increased to \/2/2, 
the lowest peak should shift upward by about 100 Hz and the highest should shift up by 
over 600 Hz. 
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CONCLUDING REMARKS 


Temperature gradients is gases affect acoustic waves by way of a spatially-variable 
sound speed. In general, the sound speed and its derivatives appear in the wave equation 
in a function which is a variable coefficient of the acoustic potential. When the gradients 
are small, the derivatives can be neglected and the function depends only on the speed of 
sound. This approximation has been called the weak gradient theory in this paper. There 
are problems where the gradients are not small, (on the scale of a typical wavelength) 
and a modified wave equation has been developed here for these cases. This modified 
equation has been called a strong gradient wave equation. The essential modification 
is the addition of the square of the sound speed gradient to the square of the circular 
frequency. The acoustic boundary conditions are altered by adding a negative imaginary 
admittance (3 T to the wall. The magnitude of this admittance is the ratio of the normal 
derivative of the sound speed to the circular frequency. 

The strong gradient theory has been applied to determine the resonant frequencies of 
rectangular cavities with temperature gradients in the direction of a one coordinate. It 
has been shown that the resonant frequencies of the cavity are changed in proportion to 
the square of the temperature gradient. Asymptotic formulas, accurate to second order in 
the temperature gradient, have been developed for the resonant frequencies and normal 
modes of the cavity. 

These formulas show that the modes with small variation in the gradient direction will 
all shift to lower resonant frequencies. Higher modes, both symmetric and antisymmetric 
in the gradient direction, may shift to lower or to higher resonant frequencies, depending 
on their zero gradient frequency. 

A conceptual experiment has been defined to study the effects of gradients on reso- 
nant frequencies. The experiment uses a two-dimensional waveguide roughly 20 cm high 
by 50 cm long to study the frequency range from 400 Hz to 4 kHz. The average tem- 
perature in the waveguide is held fixed at about the boiling point of water, 393 K while 
varying the temperature of the upper and lower surfaeo S of the waveguide. Heating the 
upper surface while cooling the lower produces a thermally-stable stratification of air with 
a nearly-constant temperature gradient. Symmetric excitation of the waveguide should 
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excite longitudinal, nearly-plane, modes whose resonant frequencies will shift downward 
with increasing temperature gradient. Antisymmetric excitation may be used to excite 
longitudinal modes whose frequencies should shift upward with increasing gradient. This 
experiment will provide a critical test of the strong gradient theory. 
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TABLE I 


Waveguide Experiment Temperatures 


L=49.14 cm, H=19.66cm, T=393.15 K 

AT/T 

0 

1/2 

y/2/2 

■s/3/2 

1 

Ti 

393.2 

294.9 

254.2 

229.9 

196.6 

t 2 

393.2 

491.4 

532.1 

563.4 

589.7 

AT,K 

0 

196.6 

278.6 

340.5 

393.2 

T , K/cm 

0 

11.80 

14.14 

17.32 

20.00 


TABLE II 


Waveguide Natural Frequencies 


L=49.69 cm, H=19.87cm, f=393.15 K 

n /l 

0 

1 

2 

3 

4 

0 

0 

1000 

2000 

3000 

4000 

1 

400 

1077 

mum 

3027 


2 

800 

WEm i 

li 

3105 


3 

1200 

1562 


3231 


4 

1600 

1887 

2561 

3400 


5 

2000 

2236 

2828 

3606 


6 

2400 

2600 

3124 

3842 


7 

2800 

2973 

3441 



8 

3200 

3353 

3774 



9 

3600 

3736 




10 

4000 

| 
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TABLE III 

Frequency Shift Constants 








TABLE V 

First Antisymmetric Mode Resonance 


fnOl > Hz 


AT/T 

0 

1/2 

mam 

■so 

1 


1000 

978 

956 

932 

909 

1 

■Bma 

1071 

1065 

1059 


2 

wmm 

1276 

1271 

1266 


3 

1562 

1563 

1564 

1565 

1566 

4 

1887 

1889 

1890 

1892 

1893 

5 

2236 

2279 

2322 

2363 

2404 

6 

2600 

2685 

2768 

2848 

2927 

7 

2973 

mmm 

3257 

3390 

3518 

8 

3353 


3784 

3982 

4170 

9 

3736 

4063 

4365 

4648 

4914 
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APPENDIX A 


Evaluation of Constants Cu 

The series for Cu converge rapidly. Since the derivatives Ajy are proportional to 
(/ — j)~ 2 , the terms in the series are proportional to (/ — j) -6 . The series for Cu can be 
arranged as the difference between a sum and a series, both of which are positive. 


C< — 
L 'll ~ 


C> = 


Cu = C> - c< 

(Al) 

Cu = 1/2 

(A2) 

2 ( t 2 + y 2 ) 2 

:+ V i 2 ' / = 3,5,... 

» (12 _ y2\ 5 

- { 3 > 

(A3) 

± ( ‘ 2+j X. < = 
^ /72 _ , 2^ 5 

=f+l,/+3,... V* J ) 

(A4) 

T (,2+J ' 2)2 i-24 

^ f/2_ 7 -2)5’ Z ’*’” 

(A5) 

Y (/ 2 + y 2 ) 2 

(A6) 


An error of about 3 5 , (less than 1%) results from taking a single term in each series, so 
that approximate values of Cu are 

C{[ = 1/2 

,2 


C| > > + i*-«y) , / = 3,5.. 

11 2/ 6 ( 2 / — l ) 5 

. (2/ 2 +2/ + l) 2 

" “ (2 1 + l ) 5 ’ ' = M,S ™ 

c > K g .. - 2 ' + < , = ,.4... 

11 ( 2 / — 1 ) 5 


(A7) 

(AS) 


C 


/ = 2,4... 

11 ( 2 / + 1) 5 


(A9) 

(A10) 

(All) 
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The difference of the positive and negative terms in the series is approximated for / > 2 by 

2/ + Q 2 _ ( 2/2 + 2/ + l) 2 _ f _ _3 _ Jl\ 

-l) 6 (2/ — l) 5 ~8 /* V 1 4/» W ’ “ 

The approximate values for Cjj can now be given as simple functions of / 




^ 193 

Cn « 

486 


M13) 


Cn ** 8/ 2 i 1 4l 2 4l A ) 1 1 2 ’ 4 ’ 6 ”- ( A14 ) 

+ / = 3 ’ 5 ’ 7 - < A15 ) 
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Figure 2. Rectangular cavity with a strong tempe 
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Figure 3. Acoustic waveguide apparatus for validating strong-gradient theory 




Figure 4. Acoustic waveguide modes 
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Figure 5. Waveguide response to symetri 
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Figure 6. Waveguide response to antisymetric excitation 
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